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ABSTRACT 

Observations and modeling for the light curve (LC) and spectra of supernova (SN) 2005bf are 
reported. This SN showed unique features: the LC had two maxima, and declined rapidly after 
the second maximum, while the spectra showed strengthening He lines whose velocity increased 
with time. The double-peaked LC can be reproduced by a double-peaked 56 Ni distribution, 
with most 56 Ni at low velocity and a small amount at high velocity. The rapid post-maximum 
decline requires a large fraction of the 7-rays to escape from the 56 Ni-dominated region, possibly 
because of low-density "holes". The presence of Balmer lines in the spectrum suggests that 
the He layer of the progenitor was substantially intact. Increasing 7-ray deposition in the He 
layer due to enhanced 7-ray escape from the 56 Ni-dominated region may explain both the delayed 
strengthening and the increasing velocity of the He lines. The SN has massive ejecta (~ 6 — 7M@), 
normal kinetic energy (~ 1.0 — 1.5 x 10 51 ergs), high peak bolometric luminosity (~ 5 x 10 42 erg 
s _1 ) for an epoch as late as ~ 40 days, and a large 56 Ni mass (~ 0.32Af Q ). These properties, 
and the presence of a small amount of H suggest that the progenitor was initially massive (M~ 
25 — 3OM ) and had lost most of its H envelope, and was possibly a WN star. The double-peaked 
56 Ni distribution suggests that the explosion may have formed jets that did not reach the He 
layer. The properties of SN 2005bf resemble those of the explosion of Cassiopeia A. 

Subject headings: supernovae: general — supernovae: individual (SN 2005bf) — stars: Wolf-Rayet — 
supernovae: individual (Cassiopeia A) 
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1. INTRODUCTION 

As the sample of well-studied supernovae (SN) 
grows, it is not surprising that rare types are ob- 
served. At first, these may seem strange and 
unique; in time they may become well-understood 
subtypes. We believe that SN 2005bf is one of 
these new types: it has unique photometric and 
spectroscopic behavior. Based on a successful ef- 
fort to model our light curve and spectra, we be- 
lieve that SN 2005bf fits into the scheme suggested 
by Nomoto et al. (1995) which places core-collapse 
SN in a sequence (HP-IIL-IIb-Ib-Ic) of increasing 
mass loss from the progenitor star. 

Discovered by Monard (2005) and Moore & Li 
(2005) on April 6, 2005 (UT) in the spiral arms 
of the SBb galaxy MCG +00-27-5, SN 2005bf was 
initially classified as a Type Ic SN (SN Ic) (Morell 
et al. 2005; Modjaz et al. 2005a). The later de- 
velopment of He lines suggested an unprecedented 
transition to Type lb (Wang & Baade 2005; Mod- 
jaz et al. 2005b). Even stranger, the light curve 
was very different from any known SN (Hamuy 
et al. 2005): a fairly rapid rise to a first peak 
was followed by a period of stalling or slow de- 
cline and by a new rise to a later, brighter peak at 
~ 40 days after explosion (Fig. 1). The bright- 
ness (Mboi ~ —18 mag) and the late epoch of 
the second peak suggest this SN ejected a large 
amount of 56 Ni. SN 2005bf does not show the 
broad lines seen in hypernovae. These properties 
make SN 2005bf a very interesting SN. 

Photometry was obtained in UBVr'i' bands at 
the 1.2m telescope at the Fred Lawrence Whip- 
ple Observatory on Mt. Hopkins (FLWO; Mod- 
jaz et al. 2005c), JHK S bands at the refur- 
bished 1.3m telescope on Mt. Hopkins (PAIRI- 
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TEL; http://pairitel.org; Modjaz et al. 2005c) 
at FLWO, BVRI bands at the Himalayan Chan- 
dra Telescope (HCT; Anupama et al. 2005), and 
BVRIJHK bands at the MAGNUM Telescope 
(Yoshii et al. 2003; Inada et al. 2005) (Fig. 1). 
Spectroscopic observations were made at FLWO, 
HCT and the SUBARU Telescope (Kawabata 
et al. 2002, 2005). We used the light curve and 
spectra to constrain models for the SN and to infer 
the properties of the progenitor star 14 . 

2. LIGHT CURVE MODELS 

The bolometric light curve (LC) was con- 
structed as in Yoshii et al. (2003; see Fig. 2). 
Synthetic bolometric LCs were computed with 
an LTE radiation hydrodynamics code and a 
gray 7-ray transfer code (Iwamoto et al. 2000). 
Electron-scattering and line opacity were consid- 
ered for radiation transport. The former used 
the ionization computed with the Saha equa- 
tion. For the latter the approximate formula 
given in Mazzali et al. (2001) was adopted, with 
modifications to account for the C, O, and He- 
rich layers. We assumed a Galactic redden- 
ing E(B — V) = 0.045, a distance modulus 
/1 = 34.5 (Schlegel, Finkbeiner, & Davis 1998; 
http://nedwww.ipac.caltech.edu/index.html; Falco 
et al. 1999), and an explosion date of 2005 March 
28 UT as inferred from the marginal detection on 
2005 March 30 UT (Moore & Li 2005). 

The LC is powered by the radioactive decay of 
56 Ni to 56 Co and 56 Fe. The theoretical LC width 
near peak depends on ejected mass M e j and ex- 
plosion kinetic energy E as M e jE~ 3 (Arnett 1982; 
Nomoto et al. 2004). The mass and distribution 
of 56 Ni are constrained by the LC brightness and 
shape, but various combinations of (M c j, E) can 
fit the LC. Spectra break this degeneracy and help 
establish the abundance distribution in the ejecta. 

The density structure used for the LC calcula- 
tion was based on the C+O star model CO138E50 
used for SN 1998bw (Nakamura et al. 2001) but 
rescaled by changing r, v, and p such that M e j 
and E are rescaled as M c j oc pr 3 and E oc pv 2 r 3 , 



After this paper was submitted and placed on astro-ph, the 
Carnegie Supernova Project and their collaborators placed 
a paper on SN 2005bf on astro-ph (Folatelli et al. 2005). It 
is beyond the scope of this Letter to compare the observa- 
tions and theoretical models in detail. 
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respectively. The He distribution is constrained 
by observational evidence and spectral fitting, as 
discussed in §3. 

The He absorption lines first appeared at a 
velocity of ~ 6,000 km s _1 . The velocity in- 
creased slowly with time, reaching ~ 7, 000 km s _1 
(Fig. 3), an unusual behavior for SN lb (see, e.g., 
Branch et al 2002). On the contrary, the velocities 
of the Fe lines declined with time. Ha, and pos- 
sibly H/3 and H7, were detected (Wang & Baade 
2005; Anupama et al. 2005), as possibly in other 
SN lb (Deng et al. 2000; Branch et al. 2002). The 
weakness of these lines suggests a small H mass 
(see §3), so we assumed that H did not affect the 
LC. Seeing traces of H suggests that the He layer 
was fairly intact before the explosion, so we in- 
cluded it in our LC computation. We assumed 
that He is mainly above v > 6,000 km s _1 , as 
suggested by the observed He line velocity The 
light curve for a model with M e j = 7Mq and 
E 5 i = _E/10 51 ergs = 2.1 is shown as a dashed 
line in Figure 2. 

The model WR star explosion of Ensman & 
Woosley (1988) produced a double-peaked LC, but 
in that case the first peak was caused by recom- 
bination of He ionized by shock heating, and is 
too narrow for the first peak of the light curve of 
SN 2005bf, which has a width of - 15 days. The 
only possibility we found to reproduce the double- 
peaked LC was to assume a distribution with most 
56 Ni at low velocity, and a small amount at high 
velocity, as follows: (X( 56 Ni), M( 56 Ni)/M Q ) = 
(0.75, 0.18) at v < 1,600 km s"\ (0.015, 0.04) 
at i) ~ 1,600 - 3,900 km s"\ (0.025, 0.07) at 
v ~ 3,900 - 6,200 km s _1 , and no 56 Ni at v > 
6, 200 km s _1 . Here X denotes the mass fraction. 
The total 56 Ni mass is M( 56 Ni) = O.29M . The 
first peak of the LC is powered by a small amount 
of high-velocity 56 Ni at v > 3,900 km s" 1 . The 
second peak is powered by 56 Ni in the low veloc- 
ity central region. The outer extent of the outer 
(v ~ 6, 200 km s" 1 ) and inner (v ~ 1, 600 km s _1 ) 
56 Ni are well determined from the time of the first 
and second peak, but the inner distribution of the 
outer 56 Ni is not well constrained, since it has little 
effect on the light curve. 

The large amount of 56 Ni (~ 0.29M©) would 
ordinarily lead to a bright 56 Co tail. The dashed 
line in Figure 2 shows a LC computed with a 7-ray 



opacity k 7 — 0.027 cm 2 g _1 (Shigeyama & Nomoto 
1990). However, the observed LC continued to 
decline rapidly, suggesting that 7-rays escape more 
easily from SN 2005bf than in the typical case. 
We used a reduced 7-ray opacity to simulate this 
situation. 

Using k 7 = 0.001cm 2 g~ 1 at v < 5,400 km s _1 
(but keeping k 7 = 0.027cm 2 g~ 1 at v > 5, 400km s _1 ), 
we obtained a set of models with (M e j/M Q , E 5 i) 
= (5, 0.6), (6, 1.0), (7, 1.3), (8, 1.7), (9, 2.3), (10, 
2.8), and (11, 3.3) that reproduce the light curve. 
The models with M j = 6 — 7M© best fit the spec- 
tra (§3). The solid line in Figure 2 shows the LC 
with M ej = 7Mq, E 51 = 1.3, M( 56 Ni) = O.32M , 
and the following 56 Ni distribution: (X( 56 Ni), 
M( 56 Ni)/M ) = (0.44, 0.22) at v < 1,600km s" 1 , 
(0.013, 0.04) at v ~ 1, 600 - 3, 900 km s" 1 , (0.03, 
0.06) at v ~ 3, 900 - 5, 400km s" 1 , and no 56 Ni at 
v > 5,400 km s _1 . 

Enhanced 7-ray escape could result from the 
presence of low-density "holes" . Such a structure 
may be produced by jets and Ray leigh- Taylor in- 
stabilities possibly caused by magnctar activity 
(see §4). While 7-rays could escape much more 
easily from the low-density holes, those 7-rays that 
are trapped in the high-density region are con- 
verted to optical photons. Thus optical photons 
are selectively produced in high density regions, 
where the optical opacity can be expected to be 
normal. 

Because of the smaller average k 7 in the 56 Ni 
region and the enhanced escape of 7-rays with re- 
spect to a model with normal k 7 , more 56 Ni is nec- 
essary to power the LC, but lower velocity 56 Ni 
(i.e., v < 5,400 km s _1 ) can reproduce the first 
peak. In other words, 56 Ni is mixed out but does 
not reach the He layer at v > 6, 000 km s _1 . Such 
a double-peaked distribution of 56 Ni might be pro- 
duced by jets that did not reach the He layer. 

3. SPECTROSCOPIC MODELS 

We computed synthetic spectra from the ex- 
plosion models described above. The comparison 
with the observed spectra is necessary to distin- 
guish among the models. We used the Monte- 
Carlo spectrum synthesis code described in Maz- 
zali & Lucy (1993), Lucy (1999), and Mazzali 
(2000). Since the code does not take into account 
non-thermal processes, which are essential to pop- 
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ulate He I levels (Lucy 1991), we introduced a 
"non-thermal factor" (Harkness et al. 1987) to 
reproduce the He lines. This parameterized fac- 
tor, denoted as /, represents the degree of depar- 
ture from the level populations computed using a 
modified nebular approximation, and is used to 
multiply all He I line opacities. 

Our synthetic spectra are compared to the ob- 
served ones in Figure 4. The bolometric luminosi- 
ties used to fit the spectra agree with the observed 
ones within 10%. The model with (M ej /M , E 51 ) 
= (7, 1.3) provides satisfactory fits for all spec- 
tra. The photosphcric velocities this model pre- 
dicts agree with those used to fit the spectra within 
20%. 

At the time of the first peak (UT April 13, 
16 days after explosion; Modjaz et al. 2005c) 
SN 2005bf exhibited SN Ic features, but actually 
both He and H can be distinguished. The feature 
near 5700A is probably He I 5876A. We obtained a 
good fit for this spectrum with a bolometric lumi- 
nosity L = 2.0 x 10 42 erg s" 1 and a photospheric 
velocity w p h = 6,200 km s _1 . The non-thermal 
factor was set to 1. The model with M e j = 7M Q 
fits better than that with M e j = 6M Q . The fea- 
ture at 6300A is reproduced as a blend of Ha and 
Si II 6355 A (Fig. 4, inset). The core of the line 
cannot be due to Si II because this would require 
a Si velocity of ~ 5,000 km s _1 , which is smaller 
than the photospheric velocity. It can, however, 
be reproduced as Ha if ~ O.O2M of H is present 
above v > 13,000 km s _1 . 

Near maximum brightness (UT May 4, 37 days 
after explosion) He lines become conspicuous. The 
model has L = 5.0 x 10 42 erg s _1 , « p h = 4,600 
km s-\ and / = 2.0 x 10 3 (only at v > 6,500 
km s _1 ). Most features, including Ca II, Fe II 
and Mg II are also well reproduced. The line near 
6300A is reproduced as Si II (6355A) at this epoch: 
the photospheric velocity is low enough for Si to 
be above the photosphere, while the H layer is too 
far above the photosphere to produce an absorp- 
tion. For a spectrum on UT May 16th (49 days, 
Kawabata et al. 2005), we also obtain a reason- 
able fit with parameters L = 3.8 x 10 42 erg s^ 1 , 
v ph = 3,800 kms" 1 , and / = 2xl0 6 (at v > 7, 200 
km s _1 ). 

These values of the "non-thermal factor" com- 
pare well with the results of detailed calculations 



(Lucy 1991), although the value at the last epoch 
appears too large. To mimic the increasing veloc- 
ity of He I A 5876, we introduced the non-thermal 
factor at increasingly high velocity at more ad- 
vanced epochs. Enhanced 7-ray escape from the 
56 Ni-dominatcd region can yield rapidly increas- 
ing 7-ray deposition in the He layer and higher He 
excitation, as confirmed by the fact that the 7-ray 
deposition rates derived from the LC calculations 
increase gradually over time at 7,200 km s _1 . The 
ratios between the deposition at 7,200 km s _1 49 
days and 16 days after the explosion is 14.5 (for 
comparison it is 4.6 at 6,200 km s _1 ). As the op- 
tical depth of the He lines increases, the region 
where they become optically thick (and where the 
non-thermal factor is applied) moves to higher ve- 
locity. 

4. CONCLUSIONS & DISCUSSION 

We have studied the properties of SN 2005bf 
by modelling the light curve and the spectra. Our 
best fit model has M ej ~ 7M and E 51 ~ 1.3. The 
ejecta consist of 56 Ni (~ O.32M ), He (~ O.4M ), 
intermediate mass elements (mainly O, Si, S), and 
a small amount of H (~ 0.02M Q ). Thus the pro- 
genitor had lost almost all its H envelope, but re- 
tained most of the He-rich layer, as in WN stars. 
The double-peaked light curve is reproduced by a 
double-peaked 56 Ni distribution: such a distribu- 
tion may be caused by jets that did not reach the 
He layer. Strong He lines were not seen in the ear- 
liest spectra because radioactive 56 Ni was too far 
from the He layer to excite He atoms. Increasing 
7-ray deposition in the He layer due to enhanced 
7-ray escape from the 56 Ni-dominated region may 
explain both the delayed strengthening and the 
increasing velocity of the He lines. 

The He core mass at explosion was Mh c = 
M oj + M cut ~ 7.5 - 8.5M . The progenitor 
was probably a WN star of main-sequence mass 
M M s ~ 25 - 3OM (Nomoto & Hashimoto 1988; 
Umeda & Nomoto 2005). The formation of a WN 
star from a star of only ~ 25M suggests that ro- 
tation may have been important (Hirschi, Meynet 
& Maeder 2005), although not sufficient to make 
SN 2005bf a hypcrnova. 

Our model requires the production of an un- 
usually large amount of 56 Ni but a normal ex- 
plosion energy. To examine whether this is pos- 
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sible, we performed hydrodynamical simulations 
of explosive nucleosynthesis for the model with 
Mms = 25M and £51 = 1.3 as in Umcda & 
Nomoto (2005). In order to produce ~ O.32M 
56 Ni, the mass cut that separates the ejecta and 
the compact remnant should be as deep as M cut ~ 
IAMq. In the spherical model no fallback oc- 
curs. This suggests that the remnant was a neu- 
tron star rather than a black-hole. The mass range 
Mms ~ 25M Q is near the transition from neu- 
tron star (SN 2005bf) to black hole formation (SN 
2002ap; Mazzali et al. 2002), the exact boundary 
depending on rotation and mass loss. 

Interestingly, the progenitor of the Cassiopeia 
A (Cas A) SN remnant was also probably a WN 
star, and its nucleosynthetic features are consis- 
tent with a ~ 25-Mq star (e.g., Fesen et al. 2005 
and references therein). The compact remnant 
of Cas A could be a magnetar (a neutron star 
with strong magnetic fields, Krause et al. 2005). 
The Cas A remnant is extremely clumpy, with 
many knots. The fastest Fe-rich 'jet' moves with 
v ~ 6, 000 km s _1 and is not mixed with N-He or 
O-rich layers (Fesen et al. 2005). These properties 
resemble our model for SN 2005bf. We might spec- 
ulate that SN 2005bf was similar to the explosion 
of Cas A and that magnetorotational effects at 
collapse and later magnetar activity produced jets 
that created extremely clumpy ejecta in the 56 Ni- 
rich layer. These jets were not energetic enough 
to reach the He layer. 

Our best-fit model with enhanced 7-ray escape 
predicts that the apparent bolometric magnitude 
of SN 2005bf as it emerges from behind the Sun in 
late 2005 will be - 20.6 (Fig. 2, inset). Without 
the enhanced 7-ray escape, this value would be 
~ 19.0 mag. SN 2005bf showed significant polar- 
ized spectral features (Wang & Baade 2005; Kawa- 
bata et al. 2005). Polarization observations (e.g., 
Kawabata et al. 2002) could be extremely useful 
to reveal the detailed distribution of the elements 
in the ejecta and the orientation of the possible 
jet, as would nebular spectra (Maeda et al. 2002; 
Mazzali et al. 2005). 
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20 40 60 80 
days since explosion (2005 March 28 UT) 

Fig. 1— (top) Observed UBVr'i'JHK s (or K) light curve of SN 2005bf. Filled symbols are data from 
the Fred Lawrence Whipple Observatory (Modjaz et al. 2005c), open symbols are data from the Himalayan 
Chandra Telescope and the MAGNUM Telescope (Anupama et al. 2005; Inada et al. 2005). (bottom) 
Observed B — V color evolution. The meaning of the symbols is the same as in the top panel. 
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-18.5 




days since explosion (2005 March 28 UT) 

Fig. 2. — The bolometric light curve constructed from FLWO (filled squares; Modjaz et al. 2005c), HCT 
(filled circles; Anupama et al. 2005), and MAGNUM (filled triangle; Inada et al. 2005) photometry. The 
contribution of near-IR light is between ~ 20% of the total at early phases to ~ 50% at late phases. Synthetic 
light curves are shown for normal (dashed) and reduced (solid) 7-ray opacities (see text). The inset shows 
the predicted LCs of each model when SN 2005bf emerges from behind the Sun in fall 2005. 
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Fig. 3. — The time evolution of the He line velocity derived from the absorption minimum of He I A 5876 
(Modjaz et al. 2005c). 
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Fig. 4.— Spectra of SN 2005bf (thick lines: 2005 April 13 - FLWO, Modjaz et al. 2005c; May 4 - HCT, 
Anupama et al. 2005; May 16 - Subaru Telescope, Kawabata et al. 2005) compared to the synthetic spectra 
(dashed lines) computed with the model (M j/M Q , -E51) — (7, 1.3). The position of He lines is shown by tick 
marks. The absorptions near 4900 and 5100A are blended with Fe II lines. The inset shows the absorption 
near 6300A in the April 13th spectrum. The model with H at v > 13,000 km s _1 (dashed line) provides 
the best fit. Thin and dotted lines show models with H in the whole ejecta and no H, respectively. 
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